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THE WARMING AND MOISTENING OF COLD AIR 
MASSES BY THE SEA 
By K. GRANT 


Summary. An empirical formula for the change in temperature of cold air crossing a 
warmer sea is based on a hypothesis adopted for the sea-to-air heat-transfer coefficient. The 
formula is tested by using routine synoptic observations made in cold airstreams reaching 
the British Isles. A formula, due to J. A. Businger, for the increase in humidity mixing 
ratio is also tested. 


Practical details are given of the use of the formulae. 


Introduction. Blackall’ has recently produced the following empirical 
formula for the increase in temperature of cold air as it crosses a warmer sea: 


T = To + (1 — exp(—12t/d)) (Ts — To). ca Oe 


(All symbols used are defined in Appendix I.) 

This paper presents the results of an attempt to find an alternative empirical 
formula for T which would give smaller root-mean-square (r.m.s.) errors than 
the above equation, and also the results of a test of Businger’s? theoretical 
formula for humidity change. 


The prognostic equations for temperature. The working hypothesis 
was adopted that the heat-transfer coefficient C, could be represented by 


cna 2) 


that is, that C, is proportional to a constant unknown power of the sea—air 
temperature difference divided by the mean air temperature. By using an 
argument similar to that used on page 66 of Blackall’s paper,’ but starting 
from the basic equation H = c,pC,, V(T; — T), this hypothesis leads to the 
prognostic equations 


T = To + [1 — exp(—pgs/d)] (Ts — To) +++ (2) 
if « = 0, and 


r= Tot [1— (14+ -5-(n- Tox) |i — 7) «(0 


ia 3 0. 
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Determination of « and f from Blackall’s synoptic data. Blackall’s 
data were examined, and, by making use of 23 of his sets of initial conditions 
covering 15 different days, revised estimates of terminal temperatures and 
dew-points were made, with a view to obtaining the temperatures right on 
the coast, and the dew-points a little way inland. It was hoped that this 
would reduce the effects of the diurnal modification of temperature over land 
and the possibly strong vertical gradient of humidity at low levels over the 
sea and coastal areas. 

For various values of « the minimum r.m.s. error and the associated optimum 
value of 8 were calculated, a range of 8 being used for each value of «. From 
these the optimum value of « and the corresponding value of f for overall 
minimum r.m.s. error in the forecast final temperature T were found. The 
sea-level air density p was obtained by assuming a constant sea-level pressure 


of 1013 mb in all cases, with T equal to (J + T>)/2. The gas constant R 
was taken as 287 J kg-? K-" (the value for dry air) and g as 9°81 ms~*. The 
results are shown in Figure 1. 

An overall minimum r.m.s. error of 0-6 degC is seen to occur at a @ 2. 
On the evidence of this analysis it was decided to adopt a value of exactly 2 
for «, and the corresponding value of 3°8 for 8 (note that the optimum value 
of 8 changes substantially for the various values of « which were tested). 

The effect of using virtual temperature was similarly explored. A minimum 
r.m.s. error of 0°57 degC at a = 2:2 was found (with B = 3:2 for « = 2). 
It was concluded that no significant benefit was to be gained by using virtual 
rather than dry-bulb temperatures, at least in these circumstances when air 
temperatures are around or below 5°C. 





(degC) 





Root-mean-square error in forecast temperature (T) 


1. 1 


2 3 


Power of (T,—T) in expression for Cy, 








FIGURE I—EFFECT ON ACCURACY OF FORECAST TEMPERATURE OF VARYING «@ IN 
Cy = B (Ts — To)*/T* 


The B which gives minimum root-mean-square error in forecast temperature is used for 
each value of «. 
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The explicit formula for temperature. When the values of « and B 
determined above are substituted, the air density p is replaced by ps/RT, the 
fetch s is expressed in terms of units of degrees of latitude (1 degree of lati- 
tude ~ 60 nautical miles ~ 111 km), and the depth d in millibars, equation (3) 
becomes 


T= To+[1— (1+ RPE) (7, — 70, me 


where N = 2-92 X 104 ps/T? (ps being expressed in millibars and Tin kelvins). 
N varies slowly with temperature and pressure as shown in Table I. 


TABLE I—VARIATION OF N WITH TEMPERATURE AND PRESSURE 


bs (mb) 
980 1000 1020 1040 


sul 


N 

1°26 

1*29 1°31 
1°32 1°34 
1°34 1°37 
1°37 1°40 
1°40 1°43 
1°43 1-46 
1°47 1-50 
1°50 2°53 
1°53 1°57 


+4++++ 


COOP NONH DOO 


Businger’s theory of humidity-mixing-ratio change. On the assump- 
tion that C, = C, and that r is constant with height through the convection 
layer, it might be expected that the formula for the evaporative flux E would 
be analogous to that for H. Businger, in Section 5.10 of Fleagle and Businger,? 
showed that this is not so, the difference being that only the temperature, not 
the mixing ratio, determines the depth of dry convection. 

Businger’s formula is 

r=1t+Z (rs — fo); -aah® 


where Z is a factor depending on the temperature increment through the 
relations 
1—Y 


Y 
and Y¥ = (T — T>)/(Ts — To). 


Z=1+ In (1 — Y) 


An interesting point here is that when Y is in the range 0-4 to o-g (as in 
practice it quite often is), Z is very nearly equal to Y — 0-2. This implies 
that if Frost’s* temperature formula T = Ty + 0-6 (Ts — To) (where To here 
is the unmodified initial surface air temperature) is accepted as being sufficiently 
accurate, then the appropriate constant in his mixing-ratio formula is not 0-6 
but 0-4. 
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Test of Businger’s formula based on the use of dew-point data. 
Businger’s formula was checked by using the dew-point and temperature data 
for the 23 occasions mentioned above. The mean error in dew-point was 
—o-36 degC and the standard deviation of the errors was 0-95 degC. These 
values compare with Blackall’s values (relative to his own terminal values) 
of —o-6 degC and 2-1 degC respectively, and show an improvement in r.m.s. 
error which is significant at the o-1 per cent level. 

During this work it was noticed that dew-points at stations actually on the 
coast often tended to be higher than those reported inland by about 1-2 degC. 
Most of this change appeared to arise close to the coast. 


An independent check of performance. Tests of the performance for 
the British Isles of equations (4) and (5), as well as that of other methods, 
were carried out by Hedge* who used independent data and trajectories 
determined from forecast charts as well as from actual charts. The results 
for actual charts are shown in Table II. 

Table II shows that, in the relatively shower-free easterlies (average sea 
fetch 180 nautical miles), the variances of the errors in the present method 
are significantly smaller than those of Frost at the 5 per cent level, but the 
reduction of Blackall’s variances by about one-third is not significant for a 
sample of only 30 trials. (If this apparent reduction in variance is true, 
however, the odds against obtaining a result which is significant at the 5 per 
cent level from 30 trials are 3 to 1.) 

In the more showery northerlies (average sea fetch 230 nautical miles), 
there is little to choose between Blackall’s method and the present method, 
though Blackall forecasts the coastal dew-points better. Frost’s method shows 
large positive mean errors but with surprisingly small standard deviations. In 
these conditions none of the methods tested really achieved a satisfactory 
performance. 

It is interesting that in Table II the mean errors in dew-point found by 
Grant’s method are —2-o0 degC for easterlies and —1-6 degC for northerlies. 


It would appear that Businger’s formula tends to forecast inland rather than 
coastal dew-points. 


Practical procedure 


(a) Forecast technique. Use the technique to forecast surface air temperatures 
and dew-points at sea or at coastal stations (also inland when vertical 
heat- and water-vapour fluxes are small) in air with near-surface temperatures 


lower than the underlying sea surface temperature, when an upwind radiosonde 
ascent is available. 


(b) Sea fetch (s). Decide on the trajectory and the representative speed V 
of the cold air. Write down the expected time of arrival of the air at your 
station (using V) and the length s of fetch over the sea in units of 1 nautical 
mile or 60 nautical miles (1 degree of latitude). 


(c) Sea temperature (Ts) and saturation mixing ratio (rs). Write down the mean 
value of these as accurately as possible. 
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(d) Depth of convection (d). On a tephigram showing the upwind station 
sounding draw the isobar of the mean sea-level pressure ps along the trajectory. 
Draw a dry adiabatic from the mean sea surface temperature and pressure 
to meet the environment curve. Draw the isobar through this meeting point. 
Write down the depth of convection in millibars (see Figure 2). 


, environment 
Je 








FIGURE 2—CONSTRUCTION FOR OBTAINING MODIFIED INITIAL SURFACE AIR 
TEMPERATURE (7) AND MODIFIED INITIAL MIXING RATIO (fo) 


A=B c<D 

DALR = Dry adiabatic lapse rate 

SALR = Saturated adiabatic lapse rate 
See Appendix I for explanation of symbols. 


(e) Modified initial air temperature (To). Draw the expected lifting conden- 

sation level (LCL) isobar at the end of the sea crossing. The best first guess 
at this is about 500 m or (fs — 60) mb. If the upper-air station is inland, 
modify the lower part of the ascent to correspond to station reports on the 
coast at the time when the air crosses it as an offshore wind. Draw the path 
followed by a parcel of air condensing at the expected LCL so that equal 
areas are enclosed between the path curve, the environment curve, and the 
isobars f; — d, as in Figure 2. The temperature at which this path curve 
intersects the surface isobars is the modified initial temperature T>. Write 
down T, and JT, — 7p. 
Note: It is necessary to modify the initial surface temperature (and mixing 
ratio) in a way such as this in order to obtain a value representative of the 
air following the vigorous mixing of the lower layers which occurs very soon 
after the air moves over the sea. 


(f) Modified initial mixing ratio (ro). Modify the dew-point curve by using 
a constant-mixing-ratio line up to the LCL and a saturated adiabatic above, 
so that the total water-vapour content between the isobars fs and ps — d is 
unchanged. This is not easily done on a tephigram since the ‘equal-areas’ 
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rule does not apply, but with experience ro should be determinable to 
within o-1 g kg~! and is given by the modified mixing ratio below the LCL 
(Figure 2). 


(g) Final temperature (T). Use the nomogram (Figure 3) to obtain the rise 
of surface air temperature AT and thence the final downwind temperature 
T= T+ AT. Figure 3 is constructed for N equal to 1-40. 


(h) Final mixing ratio (r), dew-point (Tq) and lifting condensation level. Read 
off Z to the nearest hundredth from the nomogram. The final mixing ratio 
is given by r=19+ Z(rs — 1%). Write down the corresponding surface 
dew-point 74 as determined from a tephigram. The final LCL is given 
approximately by 120 (7 — Tq) metres or 400 (T — Ta) feet. 


Notes: (1) When N (Table I) is much different from 1-40, the nomogram can 
still be used if s is multiplied by N/1-40. This elaboration will rarely be 
necessary near the British Isles. 


(2) If the sea temperature is very variable the procedure may have to be 
carried out on two or more segments serially. 


(3) Slightly greater accuracy may be achieved by successive approximation, 
using the calculated LCL, but this will rarely be worth while. 


Conclusions. The method presented in this paper for forecasting coastal 
temperatures in cold-air advection with onshore winds appears to perform 
rather better than Blackall’s method, though giving no real improvement in 
showery northerlies. Businger’s method for forecasting dew-points gives 
coastal values which are too low, often by about 2 degC, but it should be useful 
for forecasting dew-points some way inland. 


Acknowledgements. The author would like to thank Mr R. M. Blackall, 
Meteorological Office College, for providing the data, and also Dr D. J. Carson 
and Dr F. B. Smith for their useful advice, and Mr C. L. Hawson for help 
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APPENDIX I—SYMBOLS 


specific heat of dry air at constant pressure 
water-vapour transfer coefficient 

heat-transfer coefficient (= H/cppV (Ts; — T)) 

depth of convection over the sea (mb) 

upward evaporative flux of water vapour above the sea 
acceleration due to gravity 

upward flux of sensible heat above the sea 





Meteorological Magazine, 104, 1975 


S So oOo 
= | fof 





Fetch in nautical miles 
Degrees of latitude units 





N 








ao = 


mad 
On 


n" 
2 
a... 
® 
U 
“ 
® 
® 
= 
o 
® 
“oO 
.— 
“o 
en 
| 
vee 


— el 
_— ww 


4 


a) 
— 
a 
o 
— 
o 

Qa 
E 
o 
= 
a 
5 
2 
Ss 
-_ 
o 
= 
s, 
fe) 
E 
” 
> 
= 
= 
® 
— 
& 
ic) 
— 
® 

Q. 
E 
2 
fo) 
o 
wn 


NO 
Ww 





FIGURE 3—-NOMOGRAM FOR DETERMINATION OF AT From s, d, AND 7, — To 


Isopleths of A7’—— 
Isopleths of Z--- 
r=f>+ Z (re — ro) 
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ps mean sea-level pressure along trajectory 
r humidity mixing ratio of air near surface 
modified initial mixing ratio 
saturation mixing ratio at temperature T; 
gas constant for air 
sea fetch; length of completed trajectory over sea 
time of sea crossing (s/V) 
surface air temperature at sea or on coast downwind of sea 
modified initial surface air temperature 
dew-point of surface air at sea or near coast downwind of sea 
sea surface temperature 
average value of T along track 
average speed of cold air 
(T — To)/(Ts — To) . 


(r — ro) (re — 0) = 1 +A 


power of (7, — T) in expression for Cy 
dimensionless coefficient in expression for Cx 
air density near sea surface 


In (1 — Y) according to Businger® 


551.507.362.2:551.576.2 
AN IMPROVED SATELLITE NEPHANALYSIS 


By R. HARRIS and E. C. BARRETT 
(Department of Geography, University of Bristol) 


Summary. Standard satellite nephanalyses have certain limitations which suggest that 
an improvement in the present procedure would be useful. From these limitations and a list 
of the requirements of a nephanalysis procedure, an improved satellite nephanalysis on three 
layers is derived which increases the information content to between three and five times 
that of the standard nephanalysis, provides a defined set of operational guidelines for the 
analyst, and allows for a wider range of forecasting and archival usages. 


Introduction. Satellite visible imagery presents a rich array of weather 
information to the user. For many purposes a summary of this information 
is adequate, for example a snapshot summary for forecasting purposes, or 
accumulated totals and their dependent means for climatological purposes. 
Satellite nephanalyses represent a summary over space of the salient features 
of a cloud field, using a discriminant analysis-classificatory form to subdivide an 
area on the basis of its cloud characteristics (cloud type, percentage cloud 
cover, etc.). 

Although several attempts have recently been made to produce objective 
nephanalyses by computer (Decotiis and Conlan,! and Shenk and Holub?), 
there appears to have been little change in the procedure for hand-drawn 
nephanalyses since the mid sixties (Godshall*). Despite minor variations,* the 
standard nephanalysis code employed in the United Kingdom Meteorological 
Office is that operated by the former U.S. Weather Bureau, and discussed 
more fully in Barrett. This, in turn, is essentially an adaptation of the pre- 
satellite nephanalysis, defined as ‘the study of synoptic charts on which only 
clouds and weather are plotted’ (Berry et alii®), and issued to provide an input 





* For example the anticyclone symbol (U.S. Weather Bureau) = the cold pool vortex symbol 
(U.K. Meteorological Office). 
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to the forecasting services. Observations plotted on the pre-satellite neph- 
analyses were as follows: cloud type, cloud amount, precipitation, weather, 
cloud ceilings and cloud-top heights. Of these, only the first two could be 
included in satellite nephanalyses, although areas thought capable of precipitation 
are outlined as ‘synoptically significant’. Several other characteristics (e.g. size 
of features, jet-stream location, striations) have since been added. 


Limitations of standard nephanalyses. Several authors have found 
that nephanalyses provide a useful data source for cloud studies, especially in 
conventional data-sparse regions. Such studies have been both direct, e.g. 
temporal cloud summaries (Sadler*), and indirect, e.g. the relative probabilities 
and intensities of rainfall (Barrett’»*). Even towards the end of the 1960s, 
nephanalyses provided ‘the only source of global data from the weather 
satellites containing quantized cloud information’ (Clapp®). The standard 
nephanalyses do, however, contain certain inadequacies, and some of the 
more pertinent ones may be outlined as follows: 


(a) No rules governing the steps to be followed in constructing the neph- 
analysis have been established. This has resulted in a variation in 
nephanalysis procedures from country to country and from one analyst 
to another. 

(b) Only four cloud-cover categories—‘Open’, ‘Mostly Open’, ‘Mostly 
Closed’ and ‘Closed’—are recognized. Of these, two have ranges of 
20 per cent and two have ranges of 30 per cent. 

‘Synoptically significant’ cloud areas are outlined totally at the discretion 
of the analyst, without reference to any advisory guidelines. 
Cloud-type recognition is concentrated on a distinction between strati- 
form and cumuliform clouds. It is indisputable that cumulonimbus and 
cirrus are rarely identified as such in middle latitudes, whilst strato- 
cumulus is under-represented on the charts in both low and middle 
latitudes. 

There is no established minimum-size reference giving guidance as to 
which cloud patterns to include in, and which to exclude from, the 
nephanalysis. 

(f) Detail on variations of cloudiness at the ‘meso’ scale (e.g. squall-lines, 
convective-cloud patterns) is sparse. 

(g) Conventional symbols are employed even though data from a new 
source and a new view are being considered. 


Requirements of an improved nephanalysis procedure. Many of the 
criticisms that can be levelled against the standard nephanalysis arise because 
it is produced essentially as an aid to synoptic-scale forecasting. The general 
information content is consequently lower than the users of nephanalysis 
archives would like to see. In particular, it seems that any improved neph- 
analysis scheme should include more details of the cloud field at both the 
‘macro’ and ‘meso’ scales, so that it will be of more general use rather than 
simply an aid to forecasting. This is particularly important in the developing 
countries where conventional observations are sparse, and in oceanic areas 
where there are relatively few island stations. In such areas more detailed 
nephanalyses could provide much information otherwise unavailable, for 
example for reasons of cost or inaccessibility. 
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Before we propose appropriate improvements to the practice of satellite 
nephanalysis a number of objectives may be identified. These include the 
following: 


The nephanalysis information content should be as accurate and as 
complete as possible, within stated constraints of scale. 

The chart must be clear and readily understandable, irrespective of the 
degree of detail that it contains. 

The satellite-image interpretation procedure should be standardized to 
minimize operator variance. 

The scheme should be capable of completion in operational cigcum- 
stances within one hour. 

The system should be sufficiently flexible, (1) for use over areas of 
different size and (2) for application of the results to a range of secondary 
problems. 


The improved nephanalysis scheme outlined below seeks to meet these 
requirements, visible images being used initially as the data base. 


The improved nephanalysis. An increase in the amount of information 
portrayed on an improved nephanalysis raises the problem of clarity: how can 
the desired level of detail be achieved without forfeiting the simple, immediate 
visual impact of the standard nephanalysis? A first step towards solving this 
basic problem involves defining the type of analysis performed on the satellite 
imagery. This can take one of two forms: (a) a description of the cloud field, 
and (b) an interpretation of the cloud features present. If these procedures 
are separated, nephanalysis may be carried out more easily by distinguishing 
between the two methods employed. An example serves to illustrate the point. 
A roughly circular cloud feature of 6 degrees of latitude in diameter in the 
Indian Ocean, with more than go per cent of cloud cover (mainly cumulo- 
nimbus and cumulus congestus cloud types, and striations spiralling towards 
a central clear area), is a description of the appearance of a tropical storm from 
satellite altitudes. It is the interpretation of this cloud area which defines it 
as a hurricane. 

We propose that description and interpretation of cloud fields be separated 
both in stage of execution and in mode of presentation. Since description is 
the easier task to undertake, we suggest further that a three-layered neph- 
analysis, comprising two descriptive layers and one interpretative layer, would 
be even better. This readily permits an increase of some 3-5 times in the 
information content compared with the standard nephanalysis. The format 
of this improved nephanalysis comprises: 


(a) a base layer on opaque paper giving details of the geography and 
cloud cover, 
a transparent overlay carrying information on cloud type and cloud 
structure, and 
a second transparent overlay depicting an interpretation of features of 
the cloud field. 

The stages to be followed in constructing this improved nephanalysis are 
summarized in Figure 1, and an example of its operation is given in Figure 2 
(an improved nephanalysis derived from the Essa g satellite imagery for 
3 February 1970 shown in Plate I). 
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Satellite image 

















1 


Outline and denote areas of 
uniform % cloud cover in the order :- 
a) Major cloud systems 
b) Cloud areas with definite boundaries 
c) Cloud areas with indefinite boundaries 


Y 


Subdivide cloud cover areas on the 
basis of specified cloud types.: 
following the same order of analysis 


' 


Identify and portray clearly— 
defined ‘pattern and structural 
features of the cloud fields ] 
Depict recognizable synoptic 


weather systems and subsystems 
in the lower troposphere 


; 


Depict recognizable synoptic 
weather systems and subsystems 
in the upper troposphere 



































Review, and revise 
where necessary 






































\ 


Complete the analysis with additional 
descriptive terms where appropriate 


; 


Satellite nephanalysis 





























FIGURE I—STAGES IN THE CONSTRUCTION OF AN IMPROVED 
SATELLITE NEPHANALYSIS 


Description of the cloud field is the first method of analysis applied to the 
satellite photograph. The results comprise the lower two layers as follows: 


(a) The spatial distribution of cloud cover is represented by an areal 
density shading on layer 1 (Figure 2a). For this purpose, five equal 
classes of percentage overcast are used (see Figure 2d). These five 
intervals replace the previous four, thereby providing more scope for 
detailed depiction of the cloud cover. An added advantage is that each 
class possesses the same range (i.e. 20 per cent). Areas of uniform cloud 
amount are outlined in the order suggested in Figure 1, firstly definite 
and secondly indefinite systems, both bordered by solid lines. The 
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suggested minimum-size constraint operates on the basis of an area 
equivalent to a 2}-degree square, so no features larger than this should 
be omitted from the nephanalysis. All single linear features of more 
than 5 degrees of longitude or latitude in length should be included. 
Cloud type and cloud structure form the bases for the second set of 
descriptive characteristics, and are represented in symbolic form on the 
first of the transparent overlays (Figure 2b). Accurate recognition of 
cloud type is a major problem in nephanalysis so the standardization 
of recognition procedures is imperative. The scheme suggested by 
Anderson!® of identifying clouds from their brightness and texture 
characteristics provides a set of operational guidelines suitable for such 
a standardization. Where cloud type varies within an area of the same 
cloud amount, dashed lines are used to depict this change rather than 
the solid-line boundaries used for cloud cover. The remaining des- 
criptive features on the second layer are size of clouds and open spaces, 
and cloud pattern and shape. The symbols for these features are given 
in the key for the improved nephanalysis in Figure 2d. 

The interpretation of the cloud field is the final layer of the improved 
nephanalysis set (Figure 2c), and involves the identification of the 
features described in the lower two layers and other associated pheno- 
mena. In compiling the list of interpretation features shown in Figure 2d, 
information was drawn from the ESSA Technical Report entitled ‘Appli- 
cation of meteorological satellite data in analysis and forecasting’ 
(Anderson et alii"), which contains a comprehensive appraisal of the 
cloud features normally found on satellite visible imagery. One signifi- 
cant addition to the list is the identification of cloud-top height, which 
is determined from the brightness, cloud type and synoptic situation 
of a given area. The cloud-top height is then placed in one of three 
simple categories—high (above 6 km), middle (2-6 km) or low (below 
2 km) (Berry et alii®)—and portrayed appropriately (Figure 2d). 

Any combination of these three layers can be used to give the required 
level of information for a particular purpose. For example, cloud cover (layer 1) 
and cloud type (layer 2) can be used for climatological applications (especially 
in weather bureaux where the forecasting programme is predominantly 
numerical), and all three layers for synoptic forecasting (especially where more 
traditional forecasting procedures are followed). 

As the improved nephanalysis is a free-hand technique, the charts can 
be completed within one hour for ‘real-time’ operational use in synoptic 
forecasting. The optimum arrangement of the three charts at the central 
forecasting /automatic picture transmission facility is for flanking charts to be 
printed in reverse on transparent material, so that they can be folded in to 
overlay the others. For facsimile transmission to regional weather centres and 
other users the three charts could be transmitted side by side. 


Standard and improved nephanalyses: a concluding note. The im- 
proved nephanalysis introduced in this paper is an attempt to fulfil the 
requirements of a nephanalysis procedure listed above by increasing the in- 
formation content and standardizing the rules for the preparation of the 
nephanalysis. The sectionalized format of this new chart is the most significant 
contribution to the increase in the amount of information presented to the 
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FIGURE 2—-IMPROVED NEPHANALYSIS DRAWN FROM PLATE I 
(a) cloud cover (b) cloud type and structure overlay 
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LAYER 3 - INTERPRETATION OVERLAY 
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FIGURE 2—continued 
(c) interpretation overlay (d) key to the improved nephanalysis scheme 
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user. In addition, the fuller list of descriptive and interpretative features 
(Figure 2d) gives greater scope for an analysis in detail down to the ‘meso’ 
scale. The separation of (a) the method of analysis, and (b) the different types 
of information shown on each layer, gives clarity to the stages followed in an 
analysis of satellite imagery, and with this flow-chart of procedures firmly laid 
down, subjective assessment and operator variance are minimized. 

Visible imagery only is considered in the improved nephanalysis scheme 
presented here. However, an extension of the data sources to include infra-red 
imagery would be useful for increasing the information content and frequency 
of coverage provided by nephanalyses. Our work at Bristol University is 
directed now towards the development of a nephanalysis scheme for infra-red 
imagery. The biggest problem involves the diurnal temperature changes at 
the earth’s surface and in the lowest layer of the atmosphere. The biggest 
advantage stems from the twice-daily reception of the images from the satellite. 
It is to be hoped that the analytical techniques which we are designing to 
be applied by cheap objective means may be applicable either to infra-red 
or to visible satellite images. Such developments must await the future, but 
it can be affirmed now that by using infra-red and visible imagery together, 
improvements in our analysis of satellite imagery will follow. 
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551-509-324.2 
FORECASTING CONVECTIVE THUNDERSTORMS, HAIL AND 
SHOWER ACTIVITY IN THE MIDLANDS 


By N. J. ATKINS 


Summary. The parcel method, found to be widely used by forecasters because of its 
success compared with other techniques, was used in an investigation of convective activity 
in the Midlands. It is shown that, in addition to the depth of the convective layer, the amount 
of potential energy in the lower part of the layer shown on a tephigram can be used to 
determine the intensity of shower activity. An extension to the parcel method is therefore 
suggested in this paper. 


Introduction. The Saunders! test of thunderstorm forecasting showed that 
‘general practice’ gave better results than other techniques, and it was possible 
with help from several forecasters to determine what the ‘general practice’ 
was. Of the 26 forecasters who completed a questionnaire, 24 used the parcel 
method, and 10 of them used it in conjunction with another technique, which 
for seven of them included the slice method. Fourteen forecasters examined 
the size and shape of the positive area as given by the parcel method, but 
their reasons for doing this were too varied to be classified. It was decided 
therefore to analyse the relationship between weather reported and the data 
obtained from the parcel method, including the size and shape of the positive 
area shown on a tephigram. The completed questionnaire also showed that 
none of the 26 forecasters eliminated thunderstorms from a forecast because 
of the possible entrainment of dry air from the environment, but that the 
degree of such entrainment was used to judge whether convective activity 
would be widespread, scattered or isolated. During the present investigation, 
the water-vapour content of the two layers, 1000 to 700 mb and 1000 to 
850 mb were analysed, and also the dew-point separation for the same two 
layers, but no useful result was obtained. 


Period and area of the investigation. Five observing stations in the 
Midlands, namely Cardington, Wittering, Shawbury, Watnall and Elmdon 
(see Figure 1) grouped together in the Daily Weather Report were used for the 
four-year period from 1970 to 1973 inclusive. There are no radiosonde stations 
in the Midlands but normally the ascent from one or more of the nearest four 


stations (also shown in Figure 1) provided a representative temperature 
sounding. 


Shape and size of the positive area (E). A quick method was required 
to assess the shape and size of the positive area as given by the parcel method 
using the forecast maximum temperature and the forecast dew-point at the 
time of this maximum temperature. It was found that the areas could be 
represented by the sum of the separations of the parcel and environment curves 
taken at 50-mb intervals, and that it was quicker to measure them with a scale 
of millimetres than to express them in degrees read from the tephigram. 
Measurements were made on Metform 2810, and it would be necessary to 
make a conversion if a differently scaled tephigram were used. (The scale of 





Meteorological Magazine, 104, 1975 








T 
2°W 


AUGHTON 


Wittering 


Shawburye 


@Elmdon 


@ Cardington 


CAMBORNE 
| ee 


FIGURE I—-MAP SHOWING THE LOCATION OF THE FIVE REPORTING STATIONS IN 
THE MIDLANDS USED IN THE INVESTIGATION OF SHOWERS AND THUNDERSTORMS 


The positions of the nearest four radiosonde stations (4) are also shown. 




















Metform 2810 used is 10 degC = 52 mm and the energy equivalent is 
9°4 X 107% J cm~*.) The area was classified into one of three types: 


(a) with a low centre of gravity, 
(b) with a high centre of gravity, or 
(c) with a central centre of gravity. 


The positive area was divided into a lower part (£,) and a higher part (£,), 
and the classifications were made by comparing the areas. The level at which 
the area was divided was defined as the arithmetic mean of the surface pressure 
and the pressure where the parcel cut the environment curve. A measure of 
the degree of elongation of the positive area was made by linking it with the 
depth of instability at the saturated adiabatic lapse rate (SALR). Hamilton? 
has shown that the depth of instability at the SALR is a useful guide to shower 
activity. 


Weather distribution diagrams. The weather reported at five Midland 
stations was related on distribution diagrams to the depth of instability at the 
SALR, to the size of the positive area, and to the size of the partial areas 
E, and E,. These were obtained from the most representative midday 
temperature sounding. The weather reported was taken from the five stations 
as a whole so that hail at one station and a thunderstorm with rain atanother 
was Classified as a thunderstorm with hail. Occurrences of snow were included 
as rain, and different symbols were allocated to the following weather classi- 
fications: (a) nil or slight shower, (b) moderate shower, (c) heavy shower, 
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(d) moderate hail shower, (e) heavy hail shower, (f) thunderstorm with rain, 
(g) thunderstorm with hail, and (h) thunderstorm with heavy hail. It was 
necessary to include every convection day during the four-year period and 
for the purpose of this work a convection day is defined as one on which 
convection would be expected to rise above goo mb—there were 693 such 
days. All three diagrams showed an increase of shower intensity with in- 
creasing depth of instability at the SALR. There was also a trend for shower 
intensity to increase with increasing positive area, E, but this was most notice- 
able in the diagram which used E,, the lower part of the positive area. This 
showed that at a constant depth of instability the intensity of weather reported 
was to a large extent proportional to E,. Even so it was difficult to draw 
boundary lines on the diagram owing to the occurrence of several polar-air 
thunderstorms with relatively shallow instability depths. The diagram was 
therefore split into two, thereby separating the ‘cold’ days when the height 
of the 0°C isotherm was below 850 mb from the ‘warm’ days when the o°C 
isotherm was at or above 850 mb; these two diagrams are illustrated in Figures 
2 and 3. For the sake of clarity, occasions of nil or slight showers have been 
omitted from the diagrams—there were 411 such occasions. If they had been 
plotted they would mostly have occurred towards the origin of the diagram. 


Hail. It is probable that hail reached the ground in places in the Midlands 
other than the five stations used in the investigation on most days when 
thunderstorms with rain only were reported. The chance of hail falling at 
one particular place, however, increases with the lowering of the o°C isotherm 
and also with increasing E,. This is shown in Figure 4 which used all convec- 
tion days when the depth of instability at the SALR was 200 mb or more 
with the occasions of nil or slight showers again being omitted. 


Conclusions. Consideration of the models of hailstorm clouds as illustrated 
on page 159 of Ludlam,® and in Hitschfeld,* led to the working hypothesis 
that the strength of the updraught in the lower part, which must exceed a 
critical value to accelerate the hailstones upward sufficiently to carry them 
into the higher part, was likely to be more important and critical than that 
in the higher part. Hence the positive area given by the parcel method on 
a tephigram was able to provide additional information by virtue of the size 
of its lower portion and shape. 


Method. (See example in Figure 5.) 

(a) On a representative temperature sounding, modified for midday, use 
the non-superadiabatic maximum temperature and dew-point at the time 
of this maximum temperature, and draw in the parcel path from the surface 
pressure (P,) up to the point where the parcel and environment curves 
meet (P3). 

(b) Mark in the mid point of the convection layer obtained by putting 
Py = 3 (P, + Ps). 

(c) Starting at 950 mb measure in millimetres the separation between the 
parcel and environment curves at every 50-mb interval up to P,. If P, falls 
within +10 mb of a 50-mb multiple, then the measurement at the latter 
level is halved and included in the total to obtain £,. 
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FIGURE 5——-AN EXAMPLE OF THE SUGGESTED EXTENSION TO THE PARCEL METHOD 


P, (805 mb) is the arithmetic mean of the surface pressure P, (1000 mb) and the pressure 
where the path curve for the parcel cuts the environment curve at P, (610 mb). 

Ey is the sum of the measured separations at 50-mb intervals between the two curves from 

o mb up to P3. 

If P, falls within + 10 mb of a 50-mb multiple, the measurement is halved and included 
in the total. In this example (a) Ey = (16/2) + 16 + 12 + 9 = 45 mm, (b) the instability 
depth at the SALR is 930(N) — 610 = 320 mb and (c) the height of the o°C isotherm is 
895 mb. From Figures 2 and 4 of this article, thunderstorms with hail would be forecast. 


(d) Obtain depth of instability at the SALR, i.e. N — P; mb (WN is the 
condensation level). 

(e) Use values from (c) and (d) above to find the expected convection 
weather from either Figure 2 or Figure 3, depending on the height of the 
environmental o°C isotherm. An indication of the chance of hail will be 
given by £, and the height of the 0°C isotherm by using Figure 4. 


Acknowledgement. The writer wishes to thank those forecasters at 
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Cottesmore, who kindly completed the questionnaire on their method of 
forecasting thunderstorms. 
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551.589.1 (492) 2519.2 


A BIENNIAL CYCLE IN THE NUMBER OF FINE DAYS 
IN THE NETHERLANDS 


By C. J. E. SCHUURMANS 
(Koninklijk Nederlands Meteorologisch Instituut, De Bilt, The Netherlands) 


Summary. AA statistical study is made of daily weather types at De Bilt over the period 
1881-1972 which demonstrates the existence of a biennial cycle in the occurrence of warm, 
dry, and sunny weather from late spring to early autumn. 


Introduction. In 1967 Davis' published an article on the summers of 
north-west Europe. In this article he presented some good evidence of the 
existence of a biennial cycle in north-western European summers—summers 
of the odd years being generally better than summers of the even years. The 
study was mainly based on seasonal means of daily maximum temperatures 
for seven stations in the British Isles and some 30 continental stations including 
De Bilt, The Netherlands. The temperature records studied varied from 30 to 
go years, the record for De Bilt comprising 65 years. 

In Davis’s study, for each station y*-values are computed for the contingency 
tables, giving the frequencies of summers with mean daily maximum tem- 
peratures higher or lower than the preceding summer for odd and even years 
separately. In Figure 1 of the article which is reproduced here, these y?-values 
are plotted and analysed according to the significance level. From this figure 
it is concluded that the biennial cycle is best developed over the British Isles, 
southern Scandinavia and large parts of France, western Germany and Austria. 
The y*-value for De Bilt by no means reaches the 5 per cent level of significance, 
which led the author to draw an isopleth of y* = 3°84 (5 per cent level of 
significance) around the southern North Sea and its bordering areas thus 
creating an isolated region for which the biennial cycle is not in evidence. 

The present paper presents the result of a statistical study of daily weather 
types at De Bilt over the period 1881-1972. It is shown that the biennial 
cycle of summer weather is also in evidence at De Bilt. 


Weather types. Weather types have been defined in such a way that 
daily mean temperature, amount of precipitation and duration of sunshine 
determine the weather of a given day at a given place. 

The frequency distribution of daily mean temperature for a certain day 
over go years has been subdivided into terciles. The lowest tercile is called 
‘below’ (B), the middle ‘normal’ (N) and the highest ‘above’ (A). Daily 
precipitation amounts and relative duration of sunshine (per cent) have also 
been subdivided into three classes, but with fixed class limits, namely for 
precipitation: dry (D) is less than 0-3 mm, moderate (M) is 0-3-4-9 mm and 
heavy (H) is 5:0 mm or more; for sunshine: 25 per cent or less is cloudy (C), 
26-49 per cent is partial (P) and 50 per cent or more is sunny (S). 

Combinations of (B,N,A), (D,M,H) and (C,P,S) define 27 weather types. 
The frequency of occurrence of the various types varies from very frequent 
(some 15 to 20 per cent of all days in a certain season) to very rare (less than 
1 per cent). Most of the types show a rather large annual variation. 
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FIGURE I—APPROXIMATE ISOPLETHS OF I PER CENT AND 5 PER CENT LEVELS 
OF SIGNIFICANCE OF x” 


(After Davis;? see text.) 


Weather types have been determined for De Bilt for each day from 1 January 
1881 onwards.? We shall only consider here the weather type ADS—warm, 
dry and sunny weather. 


Fine days in summer. The number of fine days (ADS-days) in summer 
June, July and August) at De Bilt shows a rather large interannual variation. 
For the period 1881-1972 these numbers are given in Table I. The number of 
ADS-days varies from 1 (in 1907, 1918 and 1962) to 51 (in 1947). 

The values in Table I have been added for odd and even years separately. 
It turned out that the summers of the odd years had on the average about 4 
(or 27 per cent) more ADS-days than the summers of the even years. 


Statistical significance. Various statistical tests have been applied in 
order to investigate the statistical significance of the phenomenon. 
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TABLE I—THE NUMBER OF ADS-DAYs IN SUMMER (JUNE, JULY, AUGUST) FOR 

THE PERIOD 1881-1972 

2 8 
188- 8 10 
189- II 12 
190- 17 6 14 
1gI- 15 I 
192- 5 II 
193- 29 14 
194- 20 18 
195- 14 5 12 
196- I 15 
197- 7 


According to ‘Student’s’ t-test the difference between the mean values for 
odd and even years is significant at nearly the 1 per cent level (¢ = 2-39). 
On applying the parameter-free rank-sum test, which is more widely applicable 
than the #-test, it was also found that the hypothesis of the two samples of 
odd and even years belonging to the same distribution can be rejected at the 
3 per cent level. In the latter test the years are ranked according to the number 
of ADS-days. When equal numbers occurred, those years have been ranked 
simply in chronological order. This procedure does not affect the validity of 
the test. The rank-sum T for the odd years turned out to be 2417, which 
means that for the even years this sum (Je) must be equal to 


(N(N + 1)/2) — 2417 = 1861 (N = g2). 
The expected T (u,) and its standard deviation (o,) are given by 


N, (N, + Ne + 1) 2 N,N2 (Ny + Ne + 1) 
-= - and o,"= — 








, where in the present 


case N, = N, = 46. 

For N,, N, > 10, T is normally distributed. 

By comparison with the standard normal distribution it is found that with 
T = 2417 the area under the standard normal curve between o and 


2 = Le 


or 


= 2°17 


is 0-4850, from which we may conclude that the difference between To, (or Te) 
and p, is significant at the 1 — 2(0-4850) = 0°03, or 3 per cent level. 

Though Davis from his analysis of mean maximum temperatures could not 
find the biennial cycle to be in evidence at De Bilt, our analysis clearly shows 
that, for the indicated period, differences between summers of odd and even 
years did exist. A preliminary conclusion could furthermore be that Davis’s 
study had its limitations as to the demarcation of areas influenced by the 
biennial cycle.* 





* Davis’s method of analysis has also been criticized on statistical grounds by P. B. Wright 
(Met Mag, London, 100, 1971, pp. 301-303). 
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Fine days over the year. As will be clear from the definition, fine days 
do occur mainly in summer. In the colder seasons of course DS (dry and sunny) 
more often goes with low (B) or normal (N) temperatures. Nevertheless, the 
numbers of ADS in spring and autumn are still large enough to analyse them 
for a possible extension of the biennial cycle into these seasons. 

Table II gives for each month of the year the number of ADS-days for 
even and odd years separately. 


TABLE II—THE NUMBER OF ADS-DAys IN EACH MONTH FOR THE PERIOD 1881-1972 
Odd years Even years Difference 


January 22 26 — 4 
February 42 34 + 8 
March 114 132 —18 
April 192 + 8 
May 278 +35 
June 280 +36 
July 312 +91 
August 266 +55 
September 205 +52 
October 94 +15 
November 33 +9 
December 8 —9 


It turns out that in the autumn, especially in September, the biennial cycle 
is still influencing the number of fine days. As far as spring is concerned the 
effect seems only to be present in May. 


Synoptic-climatological considerations. It is of interest to compare 
these results with conclusions which can be drawn from the monthly patterns 
of the quasi-biennial pressure oscillation, published by Murray and Moffitt.® 
From these patterns one may infer that from May to September anomalous 
high pressure occurs in odd years, as compared with the even years, over 
western to north-western Europe. This will cause in the season concerned 
more anticyclonic conditions over The Netherlands, accompanied by rather 
light winds varying in direction between north-east to east, and south. These 
circumstances easily explain the excessive occurrence of fine days. 

In early spring and also in late autumn the pressure-difference patterns over 
Europe (still according to Murray and Moffitt) are quite different from those 
of the summer months. In general they are of such a nature that they would 
not favour more ADS-days to occur in odd years than in even years. When 
nevertheless, the figures in Table II suggest that the biennial effect is extended 
into late autumn, this could possibly be attributed to the fact that after good 
summers the North Sea is warmer than after less-good summers. 
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Analytical methods in planetary boundary-layer modelling, by R. A. Brown. 200 mm 
x 140 mm, pp. xii + 148, illus., Adam Hilger Ltd, Warner House, Folke- 
stone, Kent, 1974. Price: £8.00. 


This book cannot be recommended except perhaps as a means of gaining a 
general impression of current ideas on the origin and effect of vortex rolls in the 
boundary layer which arise from a vorticity and thermal instability of the 
quasi-Ekman layer velocity profile. These rolls occur quite frequently and 
can be important in determining the exact nature of the internal structure of the 
layer. The final third of the book is concerned with these instabilities, whilst 
the first two-thirds are devoted to a more general description of boundary- 
layer methods but with an eye very much on the needs of the final part. 

The approach is a rather mathematical one, but the value of this is largely 
lost by a profusion of errors, some typographical, others more basic. For 
example on page 47 there are four errors in the space of four equations and 
one confusing bit of notation; on page 39, the basic equation, equation (4.3), is 
apparently hopelessly incorrect, and so on. 

Even the physical concepts and arguments are often unacceptable. Section 
5-2 on the surface layer is completely misleading; it implies that it is a constant 
flux layer (the fluxes actually fall off faster there than anywhere else in the 
boundary layer) and wrongly that its depth is comparable with the depth of 
the roughness elements z, (where the basic equations are in fact quite inadequate 
and therefore cannot be applied). The description of the physical reasons for 
the existence of vortex instability in the boundary layer—potentially the most 
interesting section of the book—is no help at all to the uninitiated and no 
source of admiration for those who do understand already. 

The front cover is very attractive (showing a jet-stream cloud stretching 
across Egypt and Saudi Arabia (not Jordan!)). However, the reader has to 
search rather carefully to find those aspects relevant to the subject matter and 
the casual ‘shopper’ might be quite misled as to what this book is really about. 

F. B. SMITH 


NOTES AND NEWS 


Meteorological Office awards to captains and navigators 
of civil airlines 

Since 1954 awards have been made annually to encourage civil airline captains 
and navigators to make air reports. Suitably inscribed books are awarded to 
the captains and navigators who have provided the best series of reports 
during the year under review. Captains who have given long and meritorious 
service in the provision of air reports are considered for the award of brief-cases. 

In 1974 brief-cases were awarded to Captain D. H. F. Banton, British Airways 
(Overseas Division), and to Captain K. Mountney, British Airways (European 
Division) by the Director-General at a ceremony held in the Headquarters 
of the Meteorological Office on 2 October (see Plate IV). 
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COVOS/COMESA Conference 


A successful joint conference was held in Oxford on 24-26 September 1974 
by the French COVOS (Comité sur les Conséquences des Vols Stratosphériques) 
and the British COMESA (Committee on the Meteorological Effects of Stratos- 
pheric Aircraft) groups which have been studying the possible environmental 
effects that might arise from the operation of large numbers of supersonic 
transport aircraft in the stratosphere, and the results of their research pro- 
gramme to date were discussed. 


Co-operation in motorway fog studies 


The Meteorological Office is co-operating with the Home Office, Transport 
and Road Research Laboratory and the Atomic Energy Research Establish- 
ment in a study which is related to fog on motorways. A Meteorological Office 
transmissometer has been set up alongside the M4 motorway near Theale, 
Berks., as a reference instrument for the evaluation of a number of cheaper and 
simpler visibility devices. 


Hercules participates in GATE 


The Hercules was delivered to the Meteorological Research Flight on 3 January. 
Before its departure for Dakar, Senegal, on 24 June, in order to participate in 
the GARP Atlantic Tropical Experiment of the World Meteorological Organiz- 


ation, about two months were spent on flight tests of the aircraft, two months 
on the installation of additional experimental instrumentation, and six weeks 
on flight testing experimental installations. During the first phase of GATE, the 
aircraft flew successfully on 10 missions of 8-10 hours duration each, in addition 
to three calibration flights. In general experimental instruments have func- 
tioned satisfactorily. 


British infra-red detectors in NASA weather satellites 


What have been described by the Goddard Space Center of the National 
Aeronautics and Space Administration as the first 1-14-mile resolution pictures 
ever received from a weather satellite in synchronous orbit result from the 
development of cadmium-mercury-telluride, infra-red detector elements for use 
in a synchronous meteorological satellite (SMS) launched from Cape Canaveral 
on 15 May 1974. 

The satellite’s first work will be to watch the life cycle of short-lived storms 
that could form and die without ever being observed by satellites in a lower- 
altitude polar orbit. 

The Visible Infra-red Spin-scan Radiometer (VISSR) provides day-time 
images of nearly one quarter of the earth’s surface in 18 minutes. Pictures show 
such fine detail because each image contains 14 600 transmission lines compared 
with the normal 625-lines from standard United Kingdom television pictures. 

Maps of wind direction and force for several levels in the atmosphere are 
made from time-lapse cloud motion measurements using the infra-red channels. 
The height of the observed wind velocity is deduced from the measurement of 
cloud-top temperature. This technique has proved useful during the GATE 
experiment. 
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The detectors were produced by Mullard at Southampton. Because of the 
extreme sensitivity of the detectors to minute changes in temperature, ‘heat 
maps’ of the earth’s surface are taken by the SMS and transmitted back to 
earth to provide scientists with valuable meteorological and geological data. 
The detectors are most sensitive when they are cooled to low temperatures. 
This is achieved by a radiation cooler which allows any heat generated within 
the detector to leak into space. By use of this technique, temperatures lower 
than 100 kelvins (—173 °C) are made possible. 


Laser radar installation at Beaufort Park 


The large laser radar installation at Beaufort Park has been brought into 
operation and observations of Raman scattering from the stratosphere have 
been obtained. The results are encouraging as the number densities of air 
molecules calculated from the scattering data agree with those calculated from 
radiosonde observations in the region where the two methods of measurement 
overlap. 


Building Climatology Section set up 


A Building Climatology Section consisting of one Principal Scientific Officer 
and one Scientific Officer and funded by the Department of the Environment 
(DOE) has been set up in the Climatological Services Branch to work on 
projects specified by DOE. 


Pressure data for the British Gas Corporation 
A twice daily service of actual and forecast atmospheric pressure values has 
been supplied to the British Gas Corporation Liquid Natural Gas sites at Glen 
Mavis, Airdrie, since 27 February 1974 and Partington, Manchester, since 
1 March 1974. 


Contract for data buoy 


The Department of Trade and Industry has placed a contract for the develop- 
ment and construction of the first United Kingdom national data buoy, DB-1. 
The Office is acting as design authority for the meteorological subsystem. 


Analyses of WAMFLEX project 


Analysis of data from the WAMFLEX project has now been completed; the 
object of the project was to measure the vertical transfer of horizontal momen- 
tum over large mountains (in this case the Rockies). Wind and temperature 
data secured by the Canberra when flying at a common flight level with other 
aircraft show perfect agreement with data from those aircraft, and flux profiles 
obtained from three aircraft flying in a ‘stack’ look reasonable. Significant 
transfers of momentum flux (0-3-0-5 N/m? at most levels) occurred on three 
days. 
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Transfer of radiation work 


The transfer of all non-routine work on solar radiation and instrumentation 
from Kew Observatory to Beaufort Park was completed in January 1974. 


Computerized radiosonde equipments 


Contracts were placed with Ferranti, Wythenshawe, in March for the com- 
puterized station ground equipments and automated central calibration plant 
for the Mk 3 radiosonde system. The first ground station will be tested at 
Beaufort Park during the second half of 1975 while the central calibration plant 
will be brought into operation at Eastern Road during the same period. 


Pollution sampling over the North Sea 


In the Meteorological Office’s study of the long-range transport of industrial 
pollution, a further series of sampling flights was carried out by Meteorological 
Research Flight in dry westerlies during August and September 1973. These 
included measurements of sulphate particles as well as sulphur dioxide gas. 
The analysis of the data, which has now been completed, leads to new estimates 
of approximately one-third for the fraction of sulphur pollution deposited on 
the ground before the air crossed the east coast. In the remainder of the 
programme the opportunity was taken of making similar measurements off 
the Danish and Norwegian coasts, with the aim of obtaining information on 
the further loss by deposition over the sea. 


Computer Output Microfilm in the production of the 
Monthly Weather Report 


Commencing with the January 1974 issue, printer’s copy of the tabular material 
of the Monthly Weather Report has been in the form of Calcomp film. This new 
system will eliminate manual type-setting and proof reading. Calcomp- 
plotted charts are also being used in the production of the maps of mean air 
temperature, mean daily sunshine and 30-cm earth temperature anomaly maps 
although the isopleths are still drawn by hand. 


Study of climatic variation 


A study of climatic variation over the last few hundred years has revealed the 
inadequacies of inferring meteorological data from other associated data. For 
example a comparison of the temperature indications from the #%O content of 
a Greenland ice core and from the dates of wine harvests in northern France 
with the central England series have indicated a rather low correlation. It 
seems unlikely that reliable decadal means for the time before instrumental 
temperature data can be derived from either of these sources. 

The observed meteorological data indicate that over the last 250 years 
winter temperatures have increased slightly, but large fluctuations have 
occurred about the trend line. 
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CORRECTION 


Meteorological Magazine, September 1974, p. 266, Table II. For ‘Wokingham 
U.D.’ read ‘Woking U.D.’. 
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